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• • Abstract 

# > 

If the minimal supersymmetric standard model at scales of around 1 TeV is realized 
in nature, the total top-squark pair production cross section should be measurable at 
the CERN Large Hadron Collider. In this work we present precise predictions for this 
observable, which are based upon approximate NNLO formulas obtained using soft- 
collinear effective theory methods. 



1 Introduction 

One of the main goals of the physics program at the Large Hadron Collider (LHC) is to inves- 
tigate the existence of supersymmetric partners of the fundamental particles with masses in 
the TeV range. Within the context of the minimal supersymmetric standard model (MSSM) 
with R-parity conservation, supersymmetric particles are produced in pairs. At hadron col- 
liders, the supersymmetric particles which are expected to be most abundantly produced are 
the ones which carry color charge: squarks and gluinos. Since supersymmetry is broken, the 
mass spectrum of squarks and gluinos plays a crucial part in determining which among the 
colored supersymmetric partners is experimentally accessible. Within the context of unified 
supersymmetric theories, the scalar and gaugino masses are evolved from a common high 
scale down to the energy scale of electroweak symmetry breaking. As a consequence of large 
Yukawa and soft couplings entering the evolution of the mass parameters, the third gen- 
eration of squarks can have very different masses compared to the first two generations of 
squarks. In particular, under the assumption of "natural supersymmetry" fl], top-squarks 
are expected to be relatively light in order to reproduce the correct energy scale of elec- 
troweak symmetry breaking. Therefore, the lightest of the two supersymmetric partners of 
the top quark is usually expected to be the lightest squark in the mass spectrum. Precise 
theoretical predictions of the stop pair production cross section are instrumental in setting a 
lower bound on the lightest stop mass. Moreover, if top squarks will be discovered, accurate 
predictions of the cross section for stop pair production can be employed to determine the 
masses and other properties of these particles. For these reasons, the study of the radiative 
corrections to the production of stop pairs has already a quite long history. The calculation 
of the next-to-leading order (NLO) corrections to the cross section for stop pair production 
within the context of supersymmetric quantum chromodynamics (SUSY-QCD) was completed 
15 years ago J2l. As expected, it was found that the NLO corrections significantly decrease 
the renormalization- and factorization-scale dependences of the prediction when compared to 
the leading order (LO) calculation. Furthermore, NLO SUSY-QCD corrections increase the 
value of the cross section if the renormalization and factorization scales are chosen close to the 
value of the stop mass. The NLO SUSY-QCD corrections are implemented in the computer 
programs Prospino and Prospino2 |3|. The electroweak corrections to stop pair production 
were studied in |4l|5]. While these corrections have a quite sizable effect on the tails of the 
invariant-mass and transverse-momentum distributions, they only have a moderate impact 
on the total cross section. The emission of soft gluons accounts for a significant portion of 
the NLO SUSY-QCD corrections p], which are large. For this reason, the resummation of 
the next-to-leading logarithmic (NLL) corrections was carried out in [7|. It was found that 
these corrections increase the cross section at the LHC by up to 10% of its NLO value, while 
they further decrease the scale dependence of the prediction. In [8}jl~0] next-to-next-to-leading 
order (NNLO) threshold corrections and Coulomb corrections were derived by means of re- 
summation techniques. In these studies, the resummation is carried out in Mellin moment 
space. In the last few years, a formalism based on soft-collinear effective theory (SCET), which 
allows one to resum soft-gluon corrections directly in momentum space, was developed [IT 12 



and applied to QCD corrections for several processes of interest in collider phenomenology, 



such as Drell-Yan scattering Q3J, Higgs production 14 15 , direct photon production 16 



and recently slepton pair production 17 . A similar approach was developed independently 
in [18] , where methods of SCET and non-relativistic QCD were used to resum simultaneously 
the effects of soft and Coulomb gluons. This method was applied to squark and gluino pair 
production 18 19 , where soft and Coulomb-gluon contributions were resummed at NLL or- 
der, and to top-quark pair production 20 , where resummation up to next-to-next-to-leading 



logarithmic (NNLL) order was implemented. The soft-gluon corrections to the top-quark pair 
production cross section at the Tevatron and the LHC were also studied in 21 22 , within 
the framework developed in 
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By employing this formalism, it was possible to resum 
soft-gluon emission corrections up to NNLL accuracy, as well as to derive approximate NNLO 
formulas for the total production cross section, the pair invariant-mass distribution, the top- 
quark transverse-momentum spectrum, and the top-quark rapidity distribution. 

In fixed-order perturbation theory, differential partonic cross sections involve both singular 
distributions and regular terms. The singular distributions are functions of a "soft parameter" , 
whose precise definition depends on the kinematics. 



In 13,21 



where pair-invariant mass 
(PIM) kinematics was employed, the relevant soft parameter is a/s(1 — z), where \fs is the 
partonic center of mass energy, z = M 2 /s, and M is the pair invariant mass of the energetic 
particles produced. In this case, the soft (or partonic-threshold) region corresponds to the 
limit z — > 1. In 22 , where one-particle inclusive (1PI) kinematics was employed, the relevant 
soft parameter is S4, which is obtained by subtracting the heavy-particle mass squared from 
the invariant mass of the objects recoiling against the observed heavy particle. In this case, 
the soft gluon emission region is identified by taking the limit S4 — > 0. Detailed studies of 
Drell-Yan and top pair production processes showed that, after the convolution of the hard- 
scattering kernels with the parton luminosities, the terms which are singular in the partonic- 
threshold limit provide the numerically dominant contributions to the hadronic cross sections. 
In these processes, the singular terms typically account for more than 90% of the total NLO 
corrections. Furthermore, the dominance of the partonic-threshold regions in the calculation 
of the cross sections arises dynamically, since it is due to the strong fall-off of the parton 
luminosities outside the threshold region. This phenomenon is called "dynamical threshold 
enhancement' 
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In the partonic-threshold region, the partonic cross sections factorize into the product of 
hard functions, which contain the virtual corrections to the cross section, and soft functions, 
which account for the effects of soft gluon emissions. The singular (plus) distributions in the 
partonic cross section can be obtained by solving the renormalization group equations (RGEs) 
satisfied by the hard and soft functions. The corresponding anomalous dimensions are known 
up to NNLO 23,24 . The information extracted from the RGE and the knowledge of the 



hard and soft functions at NLO allows us to obtain approximate formulas including all of the 
singular terms up to NNLO. These approximate NNLO formulas can then be matched with 
exact NLO calculations in order to obtain precise theoretical predictions for the observables 
of interest. 

Within this context, the production of top-squark pairs in the soft-gluon emission limit 
can be studied in analogy to the production of top-quark pairs. The hard and soft functions 
are matrices in color space. The soft functions in the stop-quark and top-quark production 
processes are identical; they were evaluated in 21 22 up to NLO. In contrast, the hard 



functions for the stop-pair production process are so far unknown and must be computed with 



NLO accuracy. It is important to observe that the only SUSY parameter which appears in the 
soft functions is the mass of the produced stop quark, while all of the other SUSY parameters 
appear exclusively in the hard functions. In this work, we carry out the calculation of the hard 
functions for top-squark pair production up to NLO. By combining the NLO hard functions 
with the NLO soft functions and with the anomalous dimensions entering the RGEs satisfied 
by the various terms in the factorized cross section, it is possible to resum soft-gluon emission 
corrections up to NNLL order. Here we limit ourselves to re-expand the resummed formulas 
in order to obtain approximate NNLO formulas for the pair invariant-mass spectrum and the 
stop transverse-momentum and rapidity distributions. Although our formulas enable us to 
obtain predictions for these differential distributions, we focus our attention on the total top- 
squark production cross section, which is the observable of most immediate interest in the top 
squark searches. In fact, by integrating the approximate NNLO formulas for the differential 
distributions over the complete phase-space, we obtain predictions for the total top-squark 
pair production cross section at the LHC, and we comment on the phenomenological impact 
of the NNLO corrections arising from soft emissions. In a future work, we will evaluate the 
resummed stop pair production cross section and discuss the phenomenological impact of soft 
gluon resummation at NNLL accuracy. 



2] we introduce our notation and conventions, 
2ll|22l. Furthermore, we describe the factor- 



The paper is organized as follows: In Section 
which are very similar to the ones employed in 
ization of the stop pair production cross section in the soft limit, both in PIM kinematics and 
in 1PI kinematics. In Section |3j we discuss the calculation of the soft and hard functions up 
to NLO. In Section [4], we discuss the structure of the hard scattering kernels and we present 
approximate NNLO formulas for the stop pair production process. Predictions for the total 
top-squark pair production cross section at the LHC can be found in Section |5j together with 
an analysis of the phenomenological impact of the approximate NNLO corrections on this 
observable. Finally, we collect our conclusions in Section |6j 

2 Notation 

In this paper we study the process 

iVi(P x ) + N 2 (P 2 ) -» fifo) + i*te 4 ) + X(k) , (1) 

where N\ and N 2 indicate the incoming protons in the case of a proton-proton collider as 
the LHC, while X is an inclusive hadronic final state. In the rest of the paper, we treat 
the top squarks as on-shell particles and neglect their decay. The terms neglected in this 
approximation are of order IV /mj. At the lowest order in perturbation theory, two partonic 
channels contribute to the process in Eq. (fTl): 

g(pi) + Q(P2) -> *i(Ps) + **(P4) , 

g{pi) + g(P2) -»- tttes) + t*M , (2) 

where the momenta of the incoming partons are related to the momenta of the incoming 
hadrons by p t = XiPi (i = 1,2). The diagrams contributing to the two production channels at 
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Figure 1: Tree level diagrams contributing to the quark annihilation production channel (first 
row) and gluon fusion channel (second row). 



lowest order in QCD are shown in Figure [TJ The relevant invariants for the hadronic scattering 
are 



S = (P 1 + P 2 f , T\ = (P 1 - p 3 ) 2 - ml , U x = (Pi - pa) 2 - m\ . 
In order to describe the partonic scattering, we employ the invariants 



(3) 



s = xi% 2 S = (pi + p 2 ) 2 , h = x{T\ , U\ = x 2 Ui 



M 1 



(Ps+Pif 



s-i 



s + t x + mi = (p 4 + ky 



m? 



(4) 



In Born approximation s + 1\ + u± = 0, and consequently M 2 = s and S4 = 0. 

It is well known that the kinematics of the process allows one to define different threshold 
regions. Here we consider two different cases: PIM kinematics, in which the threshold region is 
defined by the limit s — > M 2 , and 1PI kinematics, in which the threshold region is approached 
by taking the limit S4 — > 0. Both regions were employed in the study of the top-quark pair 
production cross section and differential distributions 121 .22 .25 .26 . In particular, by working 



in PIM kinematics one can calculate the pair invariant-mass distribution, while by working 
in 1PI kinematics one can evaluate the stop transverse-momentum and rapidity distributions. 
It should be emphasized that in the PIM and 1PI threshold regions the top squarks are 
not forced to be nearly at rest, as in case of the threshold region defined by the limit j3 = 
-i/l — 4rrv[ 1 / s — > 0, which is often employed in the calculation of soft-gluon corrections to the 
total cross section (7,18, 
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In the rest of the paper, we refer to the /3 — >• limit as the 
production threshold region. 

Our goal is to employ both PIM and 1PI kinematics to obtain approximate NNLO formulas 
for the total top-squark pair production cross section. Both approaches include the numerically 
large contributions arising from the emission of soft gluons. Differences between them, and 
with the production threshold calculations, arise from the way in which different sets of power- 
suppressed corrections are treated. 



2.1 PIM kinematics 

We focus first on the case of PIM kinematics. It is convenient to introduce the quantities 



M 2 M 2 „ / 4m? 



z 



s 



S ' Ktl V * M 2 ' 



^ = V 1 --T?r- ( 5 ) 



Consequently, the PIM threshold limit s — >■ M 2 corresponds to the limit 2; — )■ 1. According to 



the QCD factorization theorem 27 , the differential cross section in M and 6 (the scattering 



angle of the top squark with respect to the beam axis in the partonic rest frame) is given by 

dW^e = JM^l d i^{l^f)c P iu, tJ ^M, C ose, H ), (6) 

i,3 Jt 

where /i/ is the factorization scale, and the sum runs over the incoming partons. In the 
following we drop the subscript PIM (and the corresponding subscript 1PI) whenever there 
is no ambiguity about the kinematic scheme employed, or when a formula applies to both 
schemes. The parton luminosities ffy are defined as the convolutions of the non-perturbative 
parton distribution functions (PDF) for the incoming partons i and j: 

f dx / v \ 

ffij (y, fJ-f) = / —fim (x, (if) fj/N 2 (-,///) = f i/Nl (y) ® fj/N 2 (y) ■ (7) 

The functions Cij in Eq. (J6J) are the hard-scattering kernels, which are related to the partonic 
cross sections and can be calculated in perturbation theory. In order not to clutter the notation, 
we do not indicate explicitly the fact that the hard-scattering kernels depend on the top-squark 
masses m^ and raf 2 , the mass trig of the first two generations of squarks and of the sbottoms 
(which we assume to be all degenerate), the top-quark mass m t , the gluino mass m g , and the 
ti-t? mixing angle a. The expansion of the C^ functions in powers of a s has the generic form 



C^ = a s 



cf + ^> + (^) ! cf + 0(a 5 



(8) 



Only the quark annihilation and gluon fusion channels contribute to C^ at lowest order in 
perturbation theory; in particular 



W-xn ,) Cf r lUl mh 



C2 = 6{l-z) 7 ^(c^-C A )l^--£ + ^), (9) 



Ml„ r \ tlUl 2m ^ 2m ' h 

(N 2 -l) \~ r t lUl A )\yM i ' M 2 t lUl 
where N = 3, and the Mandelstam invariants t\ and U\ can be written in terms of s and 9 as 

(10) 



M 2 

-—(1-/^0080), 


Ml = 

5 


M 2 



In order to calculate higher-order corrections to C qq and C gg one needs to consider virtual and 
real-emission corrections to the Born approximation. Starting at order a s , new production 
channels, such as qg — > tit\q, open up. When working in the threshold limit z — > 1, the 
calculations are simplified by the fact that there is no phase-space available for the emission 
of additional (hard) partons in the final state. Consequently, both the hard-gluon emission 
and the additional production channels are suppressed by powers of (1 — z) and can be safely 
neglected while deriving results within the partonic-threshold limit. By neglecting power- 
suppressed terms in the integrand, Eq. ^ can be rewritten as 

d 2 a 7T/% f 1 dz 



J ~z~ ^(^/)Qw(^ M 5 cos M/) 



dMd cos 9 SM 

+ ffqq (-, Vf) C m ( z > M > cos #> A*/) + ffm (-> Av) Cgg ( z > M > ~ cos e , A*/) 



'IV, 



In Eq. (11) the quark channel luminosities ff qq - and ff qq are understood to be summed over 
all light quark flavors. The two terms in the second line of Eq. ( JXIj ) differ in the fact that in 
the first term the quark (antiquark) comes from the hadron N-y (N 2 ) in Eq. (flj, while in the 
second term the quark (antiquark) comes from the hadron N 2 (Aq), respectively. The total 
cross section can be obtained by integrating over cos# in the range [— 1, 1] and over M in the 
range [2m^ , y/S] . 

In the soft-gluon emission limit z — > 1, the hard-scattering kernels Cy factor into a product 
of hard and soft functions: 

Cij(z,M, cos 9, fi f ) = Tr [H ij {M,coae,fif)S ij {y/8(l - z),M, cos 9, /i f )] + C(l - z) . (12) 

Here and in what follows we employ boldface fonts to indicate matrices in color space, such 
as the hard functions H^ and the soft functions Sij. For simplicity we drop the the top mass 
and the SUSY parameters from the list of arguments of the hard functions Hy, as well as the 
stop mass from the arguments of the soft functions S^. Throughout this paper, we work in 
the s-channel singlet-octet basis 

V C 1 ){a} ~ u aia 2 u a 3 a 4 , \^2 ) { a } ~ i a 2 a 1 i a 3 a i 5 

( r 99\ —Xa.ia.2X ( r 99\ _ • £a x a 2 c ,c (gg\ _ ja x a 2 c .c /-,o\ 

\ C 1 ){a}~ ° ° a i a i ' l C 2 ){a} ~ lJ l a 3 a 4 ' l C 3 ) { a } ~ " l a s a 4 ' \ i6 ) 

where a, represent the color index of the particle with momentum p^. We view these structures 
as basis vectors \cj) in the space of color-singlet amplitudes. Inner products in this space are 
defined through a summation over color indices as 

( Cl \ C j) = y (a)* (cj) . (14) 

{a} 

This inner product is proportional but not equal to 5jj, so the basis vectors are orthogonal 
but not orthonormal. 



A factorization formula analogous to Eq. ( 12 ) for the top-quark pair production was derived 



employing SCET and heavy-quark effective theory in 21 . A completely analogous procedure 



can be followed in order to derive Eq. (12), which is valid in the case of top-squark pair 
production. 

The hard functions are obtained from the virtual corrections and are ordinary functions 
of their arguments. The soft functions arise from the real emission of soft gluons and contain 
distributions which are singular in the z — > 1 limit. Contributions of order a™ to the soft 
functions include terms proportional to the plus distributions 



P, 



in v 



ln m (l 



m = 0,...,2n-l, (15) 



as well as terms proportional to 5(1 — z). In particular, the NLO and NNLO hard-scattering 
kernels in Eq. (Tsl) have the structure 



m=0 
3 

C £ ] = E D % P m(z) + Q^A 1 ~ *) + **(*) ■ ( 16 ) 



m=0 



i( fc ) nW n „A c( fc ) 



The functions D^j, Q Jj, and R\a depend also on cos^, M, /jj and on the heavy-particles 

masses. The coefficients D m ]ij) QqIji an d R\j can in principle be obtained from results present 
in the literature. One of the main results of this paper is the calculation of the coefficients 
D mi - (with m = 0, . . . , 3) both in the quark annihilation and gluon fusion channels. We 

can also evaluate all of the scale-dependent terms in Q jj in both channels, but due to the 
ambiguity on the choice of the normalization scale in the argument of these logarithms we 
drop part of these terms in the numerical implementation of our formulas. We will return to 
this issue below. 

2.2 1PI kinematics 

The 1PI kinematics approach allows one to describe observables in which a single parti- 
cle, rather than a pair, is considered. One can then write the top-squark rapidity (y) and 
transverse-momentum (j>t) distribution as 

~ A ^- = —^"52/ / -fi/N 1 {Xl,fJ.f)fj/N 2 {x2,fJ-f)Cipi t i j {s 4 ,S,t 1 ,Ui,fJ.f) ; (17) 

once again we will drop the subscript 1PI in most of our equations below. The expansion of 
the 1PI hard-scattering kernels Cipi in powers of a s has the same structure shown in Eq. p| 
for the PIM case. Obviously, also in this case only the qq and gg channels give non-vanishing 
contributions at lowest order in a s . The hadronic Mandelstam variables T\ and U\ are related 
to the stop rapidity and transverse momentum through the relations 



T\ = -VSm ± e- y , U x = -VSm ± e y , (18) 



where m± = > /p 2 ^ + m~ . Therefore, the variables s,s±,ti, u\, which are arguments of the 
1PI hard functions can be expressed in terms of px,y,xi, X2- The lower integration limits in 



Eq. (17) are 



X-, 



,L .-) 



xiTi 



(19) 



S + Ti l xxS + Ux 

In order to obtain the total cross section, it is necessary to integrate the double-differential 
distribution with respect to the top-squark rapidity and transverse momentum over the range 



< \y \< - In V ±f ■. 

"''"2 1-y/l-Aml/S 



<p T < 



ra 2 



(20) 



In the case of 1PI kinematics, the hard-scattering kernels in the soft-emission limit s 4 — > 
factor into a product of hard and soft functions, in analogy with Eq. ( |12| ): 

(21) 



Cij(s 4 , s', t[, u[,n) = Tr [Hij(sl, t[, u[, /j)Sij (s 4 , s', t[, u[, //)] + 0(s 4 ) . 



As emphasized in [22], the Mandelstam invariants s',t[, u[ can differ from s,t\, U\ by power 
corrections proportional to s 4 . For example, explicit results for the hard and soft functions 
can be rewritten by employing either the relation s' + t[ + u[ = or s' + t[ + u[ = s 4 . The 
difference between the two choices is due to terms suppressed by positive powers of s 4 . A 
detailed description of the way in which we deal with this ambiguity can be found in Section 4 
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As in the case of PIM kinematics, the hard and soft function are matrices in color space 
originating from virtual and soft-emission corrections, respectively. The 1PI hard functions are 
identical to the ones encountered in the case of PIM kinematics, provided that the variables 
s,tx, and Mi are written in terms of M and cos#. The soft functions are different in the PIM 
and 1PI schemes, but in both cases they are identical to the ones employed in the calculation 
of the top-quark pair production cross sections in 
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The 1PI soft functions at order a" depend on the associated plus distributions 



-Pm(<S 4 ) 

It follows that 

"ln m (s 4 /m? 



\n m (s 4 /m 2 f ) 



t\' 



s 4 



ra- 
ti 



P 

2 m 



2 

ra- 
ti 



ra 



I 

Jo 



dsA 



S.i 



9{S4) 



dsA 



ln r 



(s 4 /m|; 



s A 



[g( Si )-g(0)] + 



2n 



9(0) 



(22) 



In 



m+l 



ra 



2 

raf 



(23) 
where g is a smooth test function. At NLO and NNLO, the hard-scattering kernel in 1PI 
kinematics have the structure 



C iP = E D %Ms,) + Q { %6(s 4 ) + R%\s A ) 



m=0,l 
3 



4 2) = E d{ ZM s ±) + Qw&m + R f^) ■ 



(24) 



m=0 



As in the PIM case, the NLO coefficients D , . 



QqJj, and R\, can be in principle obtained from 
the literature. In this work we are able to derive exact expressions for the NNLO coefficients 

>(2) 



D m ij (with m = 0, . . . , 3) and the scale-dependent terms in the coefficients Qq 



(2) 



3 The hard and soft functions at NLO 

In this section we describe the calculation of the soft and hard matrices up to NLO in pertur- 
bation theory. 

3.1 Hard functions 



As it was shown in detail in 21 , the hard functions are defined in terms of the products 
of Wilson coefficients with their complex conjugates. In order to obtain the Wilson coeffi- 
cients, one matches renormalized Green's functions in SUSY-QCD with those in SCET. The 
matching procedure can be carried out by choosing arbitrary external states and infrared (IR) 
regulators. The simplest and most common procedure consists in employing on-shell states 
for the partonic processes (qq, gg) — > i\i* and in using dimensional regularization to regulate 
both IR and ultraviolet (UV) divergences. With this choice, the SCET loop diagrams vanish 
because they are scaleless. Consequently, the effective theory matrix elements are equal to 
their tree-level expressions multiplied by a UV renormalization matrix Z . The SUSY-QCD 
matrix elements are instead the virtual corrections to the scattering amplitudes of the partonic 
processes (qq,gg) — > tit*. 

Rather than implementing the matching condition for the Wilson coefficients, we directly 
move to the calculation of the the matrix elements H\^ of the hard functions projected onto a 
certain basis {|cj)} of color structures. These matrix elements, and not the Wilson coefficients, 
are needed in order to obtain the approximate NNLO formulas which represent the main result 
of the present work. Throughout this paper, we use the s-channel singlet-octet basis defined 



in Eq. (13). In order to calculate the hard-function matrix elements, we use the fact that the 



renormalized hard functions can be directly obtained from the corresponding (squared) on-shell 
QCD scattering amplitudes. The infrared poles can be removed from the QCD amplitudes by 



employing the prescription of 28 , 29 



\ M T) 



limZr 1 (e)\M ij (e))=4na s 



<E->0 



M 



ren (0) 
ij 



47T 



M 



ren (1) 



+ ... 



(25) 



where the indices ij are not summed over. Here |-My(e)) are the dimensionally regularized 
(and UV renormalized) scattering amplitudes, whose IR poles are removed by the factors 
Z^ 1 (e). The resulting finite amplitudes are expressed in terms of a s with n\ = 5 active 
flavors. The explicit form of the Z^ matrices can be obtained by employing the results 



derived in 23 24 30 



The perturbative expansion of the renormalized hard functions is defined as 



H, 



K) 



a Q 



>R 



<+° 



Arc 



< + 



(26) 



where d^ — N in the quark annihilation channel and cLr = N 2 — 1 in the gluon fusion channel. 
The desired matrix elements can then be written in terms of the renormalized QCD amplitudes 
and the color basis vectors \ci) as 



H (o)u _ I 



13 4 ( CJ | CJ ) ( CJ | CJ ) 



Cl 



M 



ren (0) 



M 



ren (0) 



i.l 



Cj 



H 



(i)i j 



4 (c/|cr) {Cj\cj) 



Cl 



M 



ren (0) 



M 



ren (1) 



Cj) + ( C/ 



M 



ren (1) 



M 



ren (0) 



Cj 



(27) 



The leading-order result for the qq channel follows from a simple calculation. In matrix 
notation, it reads 



H 



(o) 
qq 







tlMl 



1/ ti + Ux\ti 



Ml 



+ m.- 



(28 



while that for the gg channel is 



H (o) 

99 



t 



\ 



1 

N 2 

1 t\—U\ 

N M 2 


1 ti—ui 
N M 2 
(ii-wi) 2 

M 4 


1 

AT 
*1— Ml 

Af 2 


1 

N 


ti-ui 
M 2 


1 



771 J 



1t\U\ 



2M< 



m~ 



+ ^ + 4 + 2^±^ ? 



(29) 



r(l) 



In order to calculate the NLO hard matrices H^', one needs to evaluate the one-loop 
corrections to the partonic scattering amplitudes, by keeping the various color components 
separate. Although results for the corresponding one-loop diagrams interfered with the Born- 
level amplitudes are well known (2), their decomposition into the color basis is not available 
in the literature, and we therefore had to calculate it from scratch. For this purpose, we have 



used in- house routines written in the computer algebra system FORM 31 in combination with 
Reduze 



functions 34 



32,33 



The results of the calculation are expressed in terms of Passarino-Veltman 
After UV renormalization, we have derived analytic expressions for the IR 
poles, and we have evaluated the finite parts of the amplitudes numerically by employing 
the programs described in 35-37 . We have checked our results in several ways. First, by 
applying the IR renormalization factors Zij of Eq. (25), we find that the IR poles cancel 
exactly. Second, we have checked that by multiplying the one-loop hard functions with the 
corresponding tree-level soft functions and by subsequently taking the trace of the resulting 
color-space matrices, we reproduce the numerical results for the NLO virtual corrections which 
can be extracted from the code Prospino (3). 

The output of our FORM codes are too long to be presented here in explicit form. The 
explicit expressions for the NLO hard functions are coded in a Mathematica and a Fortran 
program which are included in the arXiv version of this work. In particular, these programs 
allow one to evaluate numerically the LO and NLO hard functions for an arbitrary choice of 
the input parameters. 
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3.2 Soft functions 

The soft functions are vacuum expectation values of soft Wilson-loop operators. These func- 
tions are not sensitive to the spin of the particles involved, but they depend on the color 
structure of the underlying partonic subprocesses. Consequently, the soft functions needed for 
the calculation of the top-squark pair-production cross section are precisely the same functions 
employed in the calculation of the cross section for top-quark pair production. The calculation 
of the PIM soft functions at NLO was described in detail in 



21 



while the analogous calcu- 
lation of the 1PI soft functions was carried out in 22 . For the convenience of the reader, the 



explicit results for the soft functions obtained in those two papers are collected in Appendix |A~} 

4 Structure of the hard scattering kernels 
4.1 Preliminaries 



Since the soft functions in (12) and (21) depend on plus distributions, it is more convenient 

L/2 f Sij(u,M,co8 9,fi) 



to work with the Laplace-transformed functions 

1 



Sij(L } M, cos 9, jj) 



duiexp 



uj 



e^ E jjie L 



(30) 



in PIM and similarly in the case of 1PI kinematics. It was shown in 11 13 26 that the leading 



singular terms in the hard-scattering kernels can be generated by replacing the Laplace variable 
L with a derivative d n with respect to an auxiliary variable rj, which is later set to 0. To this 
end, one defines Laplace-transformed hard-scattering kernels in PIM kinematics as 

Cij(dr!, M, cos 9, fi) — Tr[Hij(M, cos 9, n) S{ j(d n ,M, cos 9, /i)] . (31) 

For 1PI kinematics one replaces the arguments M and cos# with s', t' l5 and u[. The hard- 



scattering kernel in momentum space can then be recovered through the relations 21 , 22 26 



CpiM,ij(z,M, cos 9, n) = Cij(d v ,M, cos 9, /I) — 



Cipi,ij (s4, s', t[, u' x ,n) = Cij(d v , s', t[, u[, n) 



\HJ T(2rj) {l-z) l - 2r i 

i - 2 ^'1 1 / .s I 



7?=0 



2'/ 



r O) S4 



'm~ + SiH 



(32) 



v =o 



In order to evaluate the above formulas, one needs to employ analytic continuation to regulate 
the divergences in z — > 1 or S4 — > 0, to take the derivatives with respect to 77, and finally to 
take the limit 77 — > 0. It is possible to show that the result of this procedure can equivalently be 
obtained by implementing a series of replacement rules on the functions c^- (L, . . . ) considered 
as polynomials in a new variable L. In PIM kinematics, one must substitute 

1 -»- 5(1 - z) , 
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L -> 2P (z) + 5(1 - z)L 



M; 



L 2 -> 8Pi(z) + 4P (^)L M + 5(1 - *) L 



y iw 



-7T 



In 2 

: l-^ ; 



L 3 -> 2AP 2 (z) + 2AP 1 {z)L M + P (z) (6L 2 M - An 2 ) + 5(1 - z) (l 3 m - 2n 2 L M + 16& 



+ 



1-2 



In 2 z - 4 In ln(l - z) - 2 In zL 



M 



L 4 -» 6AP 3 (z) + 9QP 2 (z)L M + P t (z) (48L 2 7 - 32tt 2 ) + P (z) (8L 3 M - 1Qtt 2 L m + 128& 



+ 5(1 - z) [ L A M - An 2 L 2 M + 64( 3 L M + -tt 4 j + . ( _ 



27r 2 ln2-121n2ln 2 (l-z) 



12 In 2 In (1 - z)L M - 3 In zL 2 M + 6 In 2 2 In (1 - z) + 3 In 2 zL M - In 3 z 



(33) 



where L M = ln(M 2 //z 2 ), and the plus distributions P m (z) have been defined in Eq. (15). In 
1PI kinematics, one must substitute instead 

1 — > 5(s 4 ) , 

L — ► 2P (s 4 ) - 5(s 4 ) L m , 



L 2 — > 8A(S4) - 4L m P (s 4 ) + 5(s 4 ) L 



2^ 



4L 4 



3 / s 4 

I :i -> 24P 2 (s 4 ) - 24L m P!(s 4 ) + (6L 2 , - 4vr 2 ) P (s 4 ) + 5(s 4 ) (- L 3 , + 2vr 2 L m + I6C3) 
6L 4 



s.i 



(2L S - L 4 ) 



L 4 -- > 64P 3 (s 4 ) - 9QL m P 2 (s 4 ) + (A8L 2 m - 32tt 2 ) P 1 (s 4 ) + (- 8L 3 , + 16vr 2 L m + 128C 3 )p (s 4 ) 



+ 5(s 4 ) L 4 , - 4vr 2 P4 - 64C 3 ^ m + 



L A — 3L 4 L S + 3L — 27r 



15 / s 4 

,2 \ t — i„r ,2 /, ,, ,2 ,.2 



where L m = ln(/z 2 /m 2 ), L 4 = ln(l + s 4 /m 2 ), L s = ln[s|/(m~ /i 2 )], and the associated plus 
distributions P m (s 4 ) have been defined in Eq. (22). 



Since the hard and soft functions are known up to NLO, is easy to determine the NLO 
coefficient in the expansion of c in powers of a s : 

2 



ij s 



i? + £# + (£) ^ + o(c4) 



13 An l ° V4vr/ « 



(35) 



It is important to observe that the trace of the product of the LO hard function and NLO 
soft function contains the dependence of q- on L, and therefore it gives rise to the plus 
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s 


9.0 x 10 6 GeV 2 


m h 


1319.87 GeV 


tx 


-2.94979 x 10 6 GeV 2 


rriq 


1460.3 GeV 


V 


1087.17 GeV 


m t 


173.3 GeV 


m h 


1087.17 GeV 


a 


68.4° 


m~ g 


1489.98 GeV 







Table 1: Compilation of input parameters employed in the calculation of the hard-scattering 
coefficients. The angle a is the top-squark mixing angle. 

~(2) 

distributions. In order to obtain the complete coefficient q- one needs to know the hard and 
soft functions up to NNLO: 



c (2) = . - 
« d 



1 



R 



Tr 






+ Tr 



rr(0)a(2) 



Tr 



rr(2) 5 (0) 



E 

/c=0 



r (2) jk 



(36) 



The NLO hard and soft functions in this equation are know exactly. As explained in detail 
the scale-dependent part of the NNLO hard function and the L dependent part 



in 



21,22 



of the NNLO soft function can be reconstructed by exploiting the information coming from 
the RGE satisfied by these functions. To this end, one only needs the one-loop hard and soft 
matrices calculated in Section^ and in 21,22 , respectively, as well as the relevant two- and 
three-loop anomalous dimensions computed in 23 24, 38, 39 . After this information on the 



NNLO hard and soft function has been extracted, one can evaluate the coefficients c, ■■ with 



(2) 



k 



4 in Eq. ( 36 ) , as well as the scale-dependent part of Cq 



(2) 



k,ij 



4.2 Approximate NNLO formulas 

By employing the results described in the previous sections, we are able to obtain approximate 



NNLO formulas for the hard-scattering kernels in Eqs. (12) and (21). These formulas include 



the exact expressions of the coefficients multiplying the plus distributions up to NNLO, both 
in PIM and in 1PI kinematics. The complete results for these coefficients, written in terms of 
Passarino-Veltman functions, are too lengthy to be reported here. The values of the coefficients 
multiplying the various plus distributions and delta functions for arbitrary values of the input 
parameters can be extracted from the Mathematica code mentioned above, which we include 
in the arXiv submission of this work]]] As a reference for the reader, we present here the hard 
scattering kernels evaluated by setting the input parameters at the values listed in Table [Tj 
The SUSY spectrum which we chose corresponds to the benchmark point 40.2.5 in J40J; we 
employ this benchmark point for the remainder of this paper. By defining C^ = d^C^ 
(i = 0, 1, 2), in the quark annihilation channel with PIM kinematics one finds 

C { ±Jz) = 0.118673 5(1 -z), 



lr rhe Mathematica file requires the use of LoopTools 35 
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c 



(1) 

PIM,qq 



'z) = 2.53170 



.(2) 



CpiM, w -(^) = 27 - 0048 



+ 62.2067 



ln(l- 


■*Y 


1-z 


"ln 3 (l - z) 


l-z 


1 





h 1.18594 


1 


+ 

+ 


l-z_ 


+ 18.5403 


~\n 2 (l-zy 


1- 


- z 



+ 0.834825 5(1 - z) + . . . 
ln(l - z) 



56.3923 



l-z 



l-z 



(29.5324 + 82.8060) 8(1 - z) + . . . 



(37) 



where the ellipses indicate terms which are sub leading (and finite) in the z — > 1 limit. In the 
gluon fusion channel, we find instead 

4°im, s >) = 0.348572 5(1 -z), 



a 



(i) , 

PIM.ggl 



a 



(2) 
PIM,^ 1 



z) = 16.7315 


"ln(l - 
1- 


z 


+ 13.6194 

+ 


1 
l-z. 


+ 9.50848(5(1- 

+ 


z) + . . . , 


z) = 401.555 


"ln 3 (l - z) 


+ 852.324 

- + 


"ln 2 (l-z)~ 


- 389.724 

+ 


'ln(l-z)" 
l-z 


l-z 


_ l-z _ 


+ 535.481 


1 


+ 


(81.9942- 466.336) <5(1- 


-z) + .... 





(3? 



In order to completely determine the coefficients multiplying the delta functions in the NNLO 
hard-scattering kernels, one would need to know the complete NNLO hard and soft matrices. 
Since those matrices are at present not fully known at NNLO, we are only able to determine 
the scale-dependent terms in the delta-function coefficients, because those terms are governed 
by RGEs. However, since the scale-independent parts are unknown, the coefficients of the 
delta function at NNLO depend on an arbitrary second scale chosen to normalize the scale- 
dependent logarithms. In fact, for any renormalization scale \x and reference scale /j,q, one can 
always rewrite 



In 



V 


-(£ 


-Of 



(39) 



where the second term in the r.h.s. can be reabsorbed in the unknown //-independent piece. 



In Eq. (37) and Eq. (38), we chose the reference scale equal to the pair invariant mass M. 



Furthermore, in Eq. (37) and Eq. (38), we decided to separate out the contributions to the 



delta-function coefficients coming from the NNLO hard functions. In particular, in the round 
brackets multiplying the delta functions, the second number represents the contributions of 

(2) I I 

the scale-dependent terms in H\- . In our numerical analysis in Section 5 we decided to 
drop the contributions to the NNLO delta-function coefficients arising from the two-loop hard 
functions. As was observed in [22], this choice is motivated by the fact that by including these 
//-dependent terms one might induce an artificial reduction of the scale dependence, which 
can lead to an underestimated scale uncertainty. 

Similarly, in 1PI kinematics we find for the quark annihilation channel 



C im,qq^ S ^ 



0.118673 S(s 4 ) 
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a 



(i) 

lPl,qq 



'sA = 2.53170 



^lPl,qq( S 4,) 



27.0048 



+ 98.7353 




2.03883 



84.6523 



1 

«4 



+ 1.66798 S(s 4 ) + ■ 



In 2 (s 4 /m? ) 



s 4 



+ 34.0042 



ln(s 4 /m 2 ' 

s 4 



1 
s 4 



-(166.272 + 82.8060) 5(s 4 ) + 



(40) 



while for the gluon fusion channel we obtain 
C[ P \js i )= 0.348572 «5(s 4 ), 



t(i) 



C\ P \Js 4 ) = 16.7315 



A(2) 



(s 4 ) = 401.555 



+ 1511.18 




7.69240 



682.128 



1 

s 4 



+ 7.68986 <5(s 4 ) + .. 



In (s 4 /m 2 ^ 
s 4 



512.616 



ln(s 4 /m~ ' 
s 4 



1 
s 4 



(1299.25 + 466.336) <5(s 4 ) + 



(41) 



In this context, the ellipses indicate subleading terms in the s 4 — > limit. Also, in Eqs. (40) 



and (41), as in the PIM case, the scale-dependent terms are normalized to the scale s, which 



is equal to M 2 when s 4 = 0. As in the PIM case, also in Eqs. (40) and (41) the second 



number in the round brackets multiplying the delta functions in the NNLO kernels represents 
the contribution of the scale-dependent terms arising from the NNLO hard functions. 



5 Total cross section 

In this section we present a brief numerical study of the top-squark total pair production cross 
section at approximate NNLO accuracy. We stress that our approximate NNLO predictions 
include the full NLO corrections, as well as the part of the NNLO corrections arising from soft 
gluon emission, obtained by means of the procedure outlined in the previous sections. 

Following what was done in the study of the top-quark pair production total cross section, 
we employ PIM SCET and 1PI SCET kinematic schemes described in |4l]. These schemes also 
include, on top of contributions which are singular in the soft limit, NNLO terms which are 
regular in the z — > 1 (PIM SCET ) or s 4 — > (1PI SCET ) limits, and which naturally arise from 

(2) J"^ 1 J 

the SCET formalism. These terms, which are part of the functions R\a in Eqs. (|16|) and (|24|), 



do not represent the complete contribution to the NNLO cross section which is regular in the 
soft limit, since this quantity can only be obtained with a full calculation of this observable at 
NNLO accuracy. However, as shown in 21 , 22 , the regular terms appearing in the PIM SCET 
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and 1PI SCET kinematic approaches arise from the exact definition of the soft-gluon emission 
energy, and they improve the agreement between exact and approximate formulas at NLO. 
Unless otherwise stated, we present results which are obtained by averaging the ones ob- 



tained in the two kinematic schemes that we consider. In analogy with what was done in 41 
for the top-quark pair production cross section, we adopt a conservative approach and consider 
the difference between the predictions in the two kinematic schemes as an estimate of the un- 
certainty associated with the use of approximate NNLO formulas. This is justified by the fact 
that the two schemes neglect different power-suppressed terms, which are formally subleading 
but which can nevertheless have a noticeable numerical impact on the total cross section. To 
account for this uncertainty, the scale variation of the total cross section is obtained by setting 
the renormalization and factorization scales equal to each other, /ir = \xj = /i, and by varying 
this single scale between m^/2 < \i < Im^. We then look at the difference between the largest 
and smallest values obtained. In summary, the central value and perturbative uncertainties 
for the combined results at approximate NNLO accuracy are determined by employing the 
definitions 

1 / X 

Acr+ = max {cr PIM + Acr+ M , a 1Pl + Acr+J - a , 

Act m = min [a mm + Acr~ M , a 1PI + Act" J - a , (42) 

where the subscripts 1PI and PIM indicate that the corresponding quantities are evaluated in 
1PI SCET and PIM SCET kinematics, respectively, including the full set of NLO corrections and 
the contribution of the NNLO terms present in the approximate formulas for that scheme. 
As will be shown later, the total cross section is strongly dependent on the mass of the 



produced particle, m^. However, similarly to the slepton-pair production case 17 , the depen- 
dence of the total cross section on other SUSY parameters is relatively small. In order to show 
that this is indeed the case, we fix the value of the stop mass equal to m^ = 1087.17 GeV, 
corresponding to the SUSY benchmark point 40.2.5 of |40| (which is defined in Table [I]), 
and we provide predictions for the total cross section for three different sets of the remaining 
SUSY parameters. The first set of SUSY parameters coincides with the benchmark point 
40.2.5 itself. The second set is defined by choosing m t - 2 = 2639.74 GeV, rriq = 2920.61 GeV, 
nig = 2979.96 GeV, a = 136.8°. The values above are the double of the corresponding values 
of the benchmark point 40.2.5, and for this reason we label this second set of SUSY param- 
eters as "double". (The value of m t is kept fixed at 173.3 GeV.) Similarly, the third set of 
SUSY parameters (which we refer to as "half") is m t ~ 2 = 659.93 GeV, rriq = 730.15 GeV, 
rrig = 744.99 GeV, a = 34.2°. The values of the total cross sections, with the relative per- 
turbative uncertainty, for the three different SUSY parameters sets discussed above, can be 
found in Table [2} The table refers to a collider energy of 14 TeV and the PDFs employed are 



from the MSTW2008 global fit 42-44 . We observe that the numerical values for the cross 
section are very close to each other, in spite of the fact that the input SUSY parameters are 
considerably different in the three sets. This is true both at NLO and at approximate NNLO 
accuracy. 

16 



LHC 14 TeV 


MSTW2008 


SUSY point 


m h [GeV] 


1087.17 


40.2.5 


(a ± Ao- M ) NLO [pb] 


44.2±|;g x 10" 4 


"double" 


(a ± Aa /i ) NLO [pb] 


44.5+IJ x 10" 4 


"half 


(a ± Acr /i ) NLO [pb] 


42.4l 4 ;° x 10" 4 


40.2.5 


(cr ± Ao- /i ) a p prox . NNL o [pb] 


44.3+^ x 10" 4 


"double" 


(cr ± Acr /i ) a p prox . NNL o [pb] 


44.7+^j x 10" 4 


"half 


(cr ± A(T /i ) a pp r ox.NNLO [pb] 


42.3+?;^ x 10" 4 



Table 2: Stop-pair production cross sections at the LHC with \/5 = 14 TeV for three 
different sets of the SUSY parameters described in the text. The stop mass is fixed to 



m h = 1087.17 GeV. The numbers are obtained by using MSTW2008 PDFs 42-44 . 



In all of the plots and tables discussed below, the SUSY parameters are set to the values 
corresponding to the benchmark point 40.2.5 in Table [TJ This applies also to the lightest stop 
mass, with the exception of the figures in which we plot the cross section as a function of the 
stop mass or when a different choice of the stop mass is explicitly indicated. 

All of the numbers and plots are obtained by means of an in-house Fortran code, in which 
the approximate NNLO formulas are implemented. The NLO calculations, which are one of 
the elements needed to obtain predictions at approximate NNLO accuracy, have been carried 
out by modifying the public version of Prospino [3|. 

As a first step, we compare the full NLO cross section with the approximate NLO cross 
section given by the leading singular terms. To be specific, in the app roximate NLO formulas 



,(i) n(i) 



>M 



we include the coefficients -Dj J7 -, -Dq/-, and Q Qij in Eqs. (16) and (24), as well as those terms 



(i) 



ini$ 



which naturally arise in the SCET approach. The purpose of this comparison is to 
establish to what extent the leading terms in the threshold approximation reproduce the full 
cross section, or, in other words, if the dynamical threshold enhancement of the soft emission 
region takes place. This comparison is shown in Figure [2j for the case of a hadronic center of 
mass energy of 8 TeV. The two rows in the figure refer to two different choices of the PDF 
set. NLO PDFs are employed in all of the four panels. One observes that the average of the 
approximate PIM and 1PI NLO formulas reproduces very well the band obtained by varying 
the factorization scale (which is set equal to the renormalization scale) in the full NLO result. 
It is thus reasonable to expect that also the approximate NNLO formulas reproduce to a good 
extent the unknown full NNLO corrections. 

Another comparison between the scale dependence of the approximate NLO cross section 
and the full NLO cross section is shown in Figure |3j in which we plot separate curves for the 
cross section obtained in PIM SCET and 1PI SCET kinematics, and for the average of the two. 
Both for m^ = 500 GeV and for mj = 1087.17 GeV, the scale dependence of the approximate 
cross sections in the range m^/2 < fi < 1vn- tx reproduces to a good extent the scale variation 
of the full NLO cross section in the same range. 

We now turn to the discussion of the approximate NNLO results, recalling first that in 
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Figure 2: Comparison of the full NLO cross section with the one obtained at approximate 
NLO. All panels refer to the LHC with a hadronic center of mass energy of 8 TeV. In the first 
and second row we employ CT10 EBl and MSTW2008 NLO (42H44] PDFs, respectively. The 



left and right columns show different ranges in the stop mass. 

the approximate NNLO formulas we do not include the terms proportional to the Dirac delta 
functions <5(1 — z) (PIM SCET ) or <5(s 4 ) (1PI SCET ) arising from the NNLO hard functions, since 
the scale-independent parts of the NNLO hard functions are still unknown. As expected, 
the approximate NNLO predictions for the pair-production cross section show a smaller scale 
dependence than the NLO results for the same quantity. This is illustrated in Figure [4], where 
two different LHC center of mass energies (8 TeV and 14 TeV) and two different stop-quark 
masses (mf = 500 GeV and mj = 1087.17 GeV) are considered. The approximate NNLO 
line is obtained by plotting the scale dependence of the average between the PIM and 1PI 
cross sections. In order to show the effect of the approximate NNLO corrections on the scale 
dependence, both the NLO and approximate NNLO curves are plotted using MSTW2008 
NLO PDF sets. While in Figure [4] we plot exclusively the cross section scale dependence, we 
remind the reader that the perturbative uncertainty employed in all tables as well as in all 



other figures is obtained by evaluating the second and third line of Eqs. (42). Consequently, 



it reflects both the scale uncertainty and a kinematic scheme uncertainty, which is associated 
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Figure 3: Scale dependence of the PIM, 1PI, and PIM-1PI average approximate NLO cross 
sections compared to the full NLO cross section. The left panels refer to the production of a 
top squark of mass m~ tx = 500 GeV, the right panels refer to the case m t - = 1087.17 GeV. The 
two figures in the first row are obtained for LHC at v^S* = 8 TeV, the ones in the second row 
refer to \/S = 14 TeV. All of the SUSY parameters other than mi are fixed at the values of 



the benchmark point 40.2.5 40 



to the different sets of non-singular terms neglected in 1PI SCET and PIM SCET kinematics. 

A more precise assessment of the impact of the approximate NNLO corrections on the 
central value of the cross section and on the associated perturbative uncertainty can be ob- 
tained by comparing predictions for fixed values of the stop mass. This analysis is presented 
in Tables [3] and [4], [5] and [6j and [7] and [8j which refer to the LHC with center of mass ener- 
gies \/S = 7 TeV, 8 TeV, and 14 TeV, respectively. In all tables we consider two different 
values of the lightest top-squark mass: i) the value associated to the benchmark point 40.2.5, 
m^ = 1087.17 GeV, and ii) a stop mass mj = 500 GeV, close to the current lower bound 
for this particle as determined by searches at the LHC. In all tables, the first uncertainty is 
perturbative, determined as explained in Eq. (42), while the second uncertainty is obtained 



by scanning over the 90% CL sets of the corresponding PDFs and by taking into account the 
error on a s (mz)- In particular, the PDF uncertainty is determined by seeing how the average 
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Figure 4: Scale dependence of the NLO and approximate NNLO cross sections. The left 
panels refer to the production of a top squark of mass m^ = 500 GeV, the right panels refer 
to the case m^ = 1087.17 GeV. The two figures in the first row are obtained for the LHC at 
\/S = 8 TeV, the ones in the second row refer to v^S" = 14 TeV. All of the SUSY parameters 
other than m^ are set to the values of the benchmark point 40.2.5 |40|. 



of the PIM and 1PI predictions changes when evaluated with different members of the PDF 
set. The numbers for the cross sections have been obtained by employing PDFs fitted at the 
corresponding order: LO predictions are obtained by employing LO PDFs, NLO predictions 
employ NLO PDFs, and approximate NNLO predictions employ NNLO PDFs. 

In all cases listed in the tables, the inclusion of the approximate NNLO corrections reduces 
the perturbative uncertainty, when expressed as a percentage of the central value, by more 
than a factor of 2 with respect to the corresponding NLO prediction. We can summarize the 
content of the tables as follows: the scale variation in the range mj /2 < \x < 2mj can increase 
the NLO central value up to +[11, 16]% or lower it up to —[13, 18]%. At approximate NNLO, 
the scale variation can increase the central value of the cross section up to +[2, 5]% or decrease 
it up to —[5, 7]%. These considerations are valid both when one employs CT10 PDFs E5] or 
MSTW2008 PDFs. 
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LHC 7 TeV 


MSTW2008 


m h [GcV] 


500 


1087.17 


(a ± A(i M ± A pdf ) LO [pb] 


34 4+15.8+3.8 iq-3 
°^-^-10.0-3.6 A ±u 


38.8+^™ x 10 - 6 


(a ± Aa M ± A pdf+Q J NLO [pb] 


46.2^-1+6-6 x 10 -3 


49.6+|?^ x 10- 6 


(c ± Acr^j ± A pd f +Qs ) appi . ox . nnlo [pb] 


46.2+™;° x 10- 3 


52.8-tlitlH x 10- 6 


-K'nLO 


1.34 


1.28 


-f^approx. NNLO 


1.34 


1.36 



Table 3: Stop-pair production cross section for two different values of mj at the LHC with 
v5 = 7 TeV. The numbers are obtained by using MSTW2008 PDFs. Here and in the following 
tables, all of the SUSY parameters (with the exception of TOf ) are fixed at the values prescribed 
by the benchmark point 40.2.5 |40). 



LHC 7 TeV 


CT10 


m h [GeV] 


500 


1087.17 


(a ± Aff M ± A pdf+a J L0 [pb] 


30. 1+™ 1 x 10- 3 


36.7^I;|+?g;I x 10" 6 


(a ± Aa M ± Apdf+ajNLO [pb] 


45.3i^^ 1 i° x 10- 3 


58.4 + ? 3 t1 2 9 4 x 10- 6 


(c ± Acr M ± A pd f +Q J approx . nnlo [pb] 


46.7±£Jt& 8 x 10- 3 


5i.7i 3 :S 6 x io- 6 


-K'nLO 


1.50 


1.59 


-f^approx. NNLO 


1.55 


1.41 



Table 4: Stop-pair production cross section for two different values of mj at the LHC with 
iS = 7 TeV. The numbers are obtained by using CT10 PDFs. 



In almost all cases illustrated in the tables, the PDF and a s uncertainty grows marginally in 
the approximate NNLO predictions with respect to the NLO predictions. Another way to look 
at the PDF and a s uncertainty is shown in Figure [5j where this uncertainty band is plotted 
as function of the top squark mass in the range m% 1 G [200, 2000] GeV at the LHC with center 
of mass energy of 8 TeV. The left panel refers to the case in which CT10 PDFs are employed, 
while the right panel refers to the case in which the PDFs employed are MSTW2008. One sees 
that both bands become larger for large stop masses. The approximate NNLO band in the 
left panel is almost everywhere inside the NLO band, while in the right panel the approximate 
NNLO band is larger than the NLO band. However, the bands obtained by using CT10 PDFs 
remain larger than the ones obtained when using MSTW2008 PDFs. 

Tables [3] to [8] also include the values for the K factors at NLO and approximate NNLO, 
which are both normalized to the LO cross section, i.e. 



K 



NLO 



CNLO 
0"LO 



K.. 



approx. NNLO 



°"approx.NNLO 
0"LO 



(43) 



The NLO K factors tend to be slightly larger when CT10 PDFs rather than MSTW2008 
PDFs are employed (roughly 1.5 vs. 1.3), but they are not very sensitive to the collider center 
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LHC 8 TeV 


MSTW2008 


m h [GcV] 


500 


1087.17 


(<7 ± Aa M ± A pdf ) LO [pb] 


fi i 7 +27.3+6.1 iq-3 
0i -' -17.5-6.0 X J-U 


lLBttitli x lO" 5 


(o- ± Aa M ± A pdf+a J NLO [pb] 


83.4t^;^° 8 6 x lO" 3 


U.7 + _lt-lI x lO" 5 


(°" =t Aa /t ± A pd f +Q , s ) a p prox . nnlo [pb] 


83.21; 3 ;™ x 10- 3 


15.3i?; 3 t^ x lO" 5 


-K'nLO 


1.35 


1.29 


-f^approx. NNLO 


1.35 


1.34 



Table 5: Stop-pair production cross section for two different values of mj at the LHC with 
f<3 = 8 TeV. The numbers are obtained by using MSTW2008 PDFs. 



LHC 8 TeV 


CT10 


m h [GeV] 


500 


1087.17 


(a ± Aa M ± A pdf+a J LO [pb] 


54.0^;™ x 10- 3 


10.6l|-SlI1 x IO" 5 


(a ± Acr^ ± Apdf+ajNLO [pb] 


so.g+^t 1 ^ 6 ! x io- 3 


16.5ii 3 ^°3 4 x IO- 5 


(C ± AfT^ ± Apdf+ojapprox. NNLO [pb] 


83.6i 3 |^^ x 10- 3 


15.2+°; 3 +^ x IO- 5 


-K'nLO 


1.50 


1.56 


-f^approx. NNLO 


1.55 


1.44 



Table 6: Stop-pair production cross section for two different values of m^ at the LHC with 
'S = 8 TeV. The numbers are obtained by using CT10 PDFs. 



of mass energy or to the mass of the top squark. The ratio i^ apP rox. nnlo/^nlo ranges from 
0.88 to 1.06, therefore the approximate NNLO corrections have only a moderate impact on 
the central value of the NLO cross section]^] For top-squark masses smaller than ~ 1 TeV, 
the central value for the approximate NNLO cross section falls well within the NLO scale 
uncertainty band. 

Finally, Figure [6] shows the cross section as a function of the top-squark mass up to m^ = 
2 TeV for the LHC at 7, 8 and 14 TeV center of mass energy. In all cases, the bands represent 
the residual perturbative uncertainties, obtained as explained above. In Figure [6] we employ 
CT10 PDFs. One sees that, for large values of the stop mass, the approximate NNLO band is 
below the NLO band at 7 and 8 TeV center of mass energy, while, for the LHC at 14 TeV, the 
approximate NNLO band overlaps with the lower part of the NLO band. For comparison, in 
Figure [7] we repeat the same analysis by employing MSTW2008 PDFs for the LHC at 8 TeV. 
In this case, for large stop masses the approximate NNLO scale uncertainty bands tend to be 
slightly above the NLO bands. 

We conclude this Section by comparing our results with the ones presented in J7], which 
have NLO+NLL accuracy. In Tables M and 10 we show the results obtained for the input 



2 In this numerical analysis we use NNLO PDFs together with our approximate NNLO results for the hard- 
scattering kernels. One could also make a different choice and use NLO PDFs with the approximate NNLO 
formulas, in that case the impact on the central values of our predictions would be bigger. 
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LHC 14 TeV 


MSTW2008 


m h [GeV] 


500 


1087.17 


(<7 ± Aa M ± A pdf ) LO [pb] 


48.3±i|-t±ii x 10- 2 


33.5™4 3 6 7 x 10" 4 


(o- ± Aa M ± A pdf+a J NLO [pb] 


66.4i™; 2 x 10- 2 


44.2±|;g±|;f x 10" 4 


(°" =t Aa /t ± A pd f +Q , s ) a p prox . nnlo [pb] 


65.7±1;|±|;| x 10" 2 


U.StlitH x 10" 4 


-K'nLO 


1.38 


1.32 


-f^approx. NNLO 


1.36 


1.32 



Table 7: Stop-pair production cross section for two different values of mj at the LHC with 
^S = 14 TeV. The numbers are obtained by using MSTW2008 PDFs. 



LHC 14 TeV 


CT10 


m il [GeV] 


500 


1087.17 


(a ± Acr M ± A pdf+Q J LO [pb] 


42.6+^+fs x lO" 2 


30.1+^ 3 4 7 2 8 x lO" 4 


(er ± A<t m ± A pdf+Qg ) NLO [pb] 


63.2±™+™ x 10- 2 


44.1+^ 1 1 7 x 10- 4 


(c ± Act^ ± Ap d f +Q , s ) a p prox . nnlo [pb] 


65.9±m 2 x 10- 2 


44.6+^J:^ 1 x lO" 4 


-^NLO 


1.48 


1.47 


-fdpprox. NNLO 


1.55 


1.48 



Table 8: Stop-pair production cross section for two different values of mi at the LHC with 
'S — 14 TeV. The numbers are obtained by using CT10 PDFs. 



parameters employed in [7], which coincide with the SPSla' benchmark scenario of 47 , as 
well as for two different values of the stop mass. Table M refers to the LHC with 7 TeV center 



of mass energy while Table 10 refers to the LHC with a center of mass energy of 14 TeV. 
We checked that, as expected, the NLO results in [7J coincide with the ones obtained with 
our modified version of Prospino. The central values of our approximate NNLO predictions 
for the total cross section are in very good agreement with the ones obtained by means of 
NLO+NLL formulas in [7|. However, the perturbative uncertainties at approximate NNLO 
accuracy appear to be smaller that the NLO+NLL scale uncertainties. 



6 Conclusions 

Supersymmetry is certainly one of the best motivated scenarios for physics beyond the stan- 
dard model. If supersymmetry is broken just above the electroweak scale, the supersymmetric 
partners of the known particles are expected to have masses of the order of 1 TeV, and they 
could soon be observed at the LHC. For these reasons it is important to have precise theoretical 
predictions for the production cross sections of supersymmetric particles at hadron colliders. 
Top squarks are expected to be among the lightest colored supersymmetric particles, and 
consequently they might be the supersymmetric partners which are most easily accessible at 
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Figure 5: PDF and a s uncertainties on the total stop production cross section at the LHC 
with VS = 8 TeV. 



CTEQ6.6 PDFs 



{a ± Ao-jJnlq+nll [P b l from 



(a± Act 



M/approx 



. NNLO [pb], this work 



m il = 100 GeV 



390+^ 



oqo+39 



m ~ = 400 GeV 



209+ ' 018 



0.215" 



-0.010 
-0.012 



Table 9: Comparison between the NLO+NLL results of |7| a nd the approximate NNLO results 
of the present work. The table refers to the LHC with \^S = 7 TeV; the PDFs employed are 
from CTEQ6.6 46 . We report only the perturbative uncertainties. 



the LHC. In this paper we have employed effective field theory techniques in order to improve 
existing NLO calculations for stop-pair production by evaluating higher-order perturbative 
corrections arising from soft-gluon emissions. 

In particular, we have adopted a framework which allows one to resum large logarithmic 
corrections at NNLL accuracy in the stop pair production process. The resummation can 
be carried out in two different kinematic schemes, referred to as PIM and 1PI. In principle, 
this fact enables one to obtain predictions for both the pair invariant mass distribution and 
the top-squark px and rapidity distributions, as well as for the total cross section, which is 
at the moment the observable of primary interest in phenomenological studies. Furthermore, 
by re-expanding the NNLL resummation formulas, it is possible to obtain predictions for the 
cross section which have approximate NNLO accuracy in fixed-order perturbation theory. The 
evaluation of these approximate NNLO formulas represents the main goal of this work, where 
we obtained analytic expressions for all of the coefficients multiplying the plus distributions 
in the variables (1 — z) and S4 in the NNLO hard-scattering kernels for the double differential 
distributions in the two kinematic schemes considered. This was made possible by a complete 
computation of the hard-function matrices in both production channels at NLO accuracy, the 
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Figure 6: Mass scans with CTIO PDFs for the LHC with v^ = 7, 8, and 14 TeV. The bands 
represent the perturbative scale uncertainties at NLO and NNLO. 
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Figure 7: Mass scans with MSTW2008 PDFs for the LHC at VS = 8 TeV. The bands represent 
the perturbative scale uncertainties at NLO and NNLO. 



CTEQ6.6 PDFs 


m h = 100 GeV 


m k = 400 GeV 


(a ± Act m ) nlo+nll [pb] from 


7 




L65 - :i6 x 1Q3 


o i q+0.20 
z ' ly -0.19 


(a ± A<7 At ) approx . NNL0 [pb], this work 


1.64^ x 10 3 


9 99+O.I3 
z - zz -0.10 



Table 10: Comparison between the NLO+NLL results of [7] and the approximate NNLO 
results of the present work. The table refers to the LHC with V5 = 14 TeV; the PDFs 
employed are from CTEQ6.6. We report only the perturbative uncertainties. 

calculation of which represents the main technical result of the present work. 

The impact of the approximate NNLO corrections on the stop pair production cross section 
has been examined. In order to obtain the best possible predictions and a better control on the 
neglected subleading terms, we have averaged the approximate NNLO results obtained in PIM 
and 1PI kinematics and matched them with exact fixed-order NLO calculations. As in the 
case of slepton-pair production studied in 17 , we found that the total cross section depends 



strongly on the mass of the produced particles, while the dependence on the remaining SUSY 
parameters, such as the masses of other supersymmetric particles, is rather weak. We have 
found that including the approximate NNLO predictions for the pair production cross section 
reduces the perturbative uncertainty by more than a factor of two with respect to the NLO 
results. We stress that we explicitly accounted for the kinematic scheme uncertainty arising 
from the use of approximate NNLO formulas and we combined it with the scale uncertainty, 
as explained in Section [5j It was found that the approximate NNLO corrections have only 
a moderate impact on the central value of the NLO total cross section. We conclude by 
observing that, after the NNLO corrections are included, the main theoretical uncertainty on 
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the total cross section comes from the PDFs. This is particular evident when large values 
of the stop mass are considered. Therefore, the results presented in this paper allow one to 
improve the precision of the stop-pair total cross section predictions by reducing their scale 
dependence to the extent that it becomes considerably smaller than PDF uncertainty. 
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A NLO soft functions 

In this appendix we collect the results of the calculation of the NLO soft function in PIM 



kinematics 21 and in 1PI kinematics 22 



It was proven in 13 that the Laplace-transformed soft function defined in Eq. (30) is 



related to the position-space Wilson loop 48 . In particular, the soft function in position 
space is defined as (see Eq. (39) in [2~T] ) 

W{x t ») = ^(0|T[O s t(a;)]T[O s (0)]|0), (44) 

where T (T) indicates time ordering (anti-time ordering), while O s is the operator 

O s (x) = [S n S n S m S V4 \(x) . (45) 

The four soft Wilson lines S are oriented along the directions of the incoming partons (n 
and n), and along the four- velocities of the outgoing top squarks (t>3 and v±). The precise 
definition of the soft Wilson lines can be found in Eqs. (22) in [21]. The Laplace transformed 
soft functions can be obtained from the soft functions in position space through the relation 

HL^) = W^x = ^ rj - 2 ^y (46) 

Here and below, we have omitted the dependence of the soft functions on the PIM or 1PI 
kinematic variables and on the heavy particle masses, as well as the subscripts qq or gg 
indicating the channel. The expansion of s in powers of the strong coupling constant is 

s = s <o> + ^ s a> + (£|)> + 0(a 3 ) . (47) 
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At leading order the soft functions are the same both in PIM and 1PI kinematics: 
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(48 



The bare soft function at NLO in position space can be written as 

w£J(e,s , A*) = ]T 1^(6, x ,/i) , (k = PIM, 1PI) . 



(49) 



The matrices Wij are related to the products of color generators and are the same for both 
kinematics. In the quark annihilation channel they are 



U7 12 — "^34 



C F (AN 2 
'an I -1 



w 



n 
:« 



w 



ii 

44 



C F 2N 
~2~ I C J ' 



It; 



13 



it; 



24 



C F 1 

'~2~ll 2C F - W J 



it; 



14 



it; 



ii 

23 



C F / -TV' 
2AM -TV 1 



while for the gluon fusion channel one finds 
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(50) 



toff = toff 



/ -AN \ 

_4AT iV 2 -(iV 2 -4) 

\ -(iV 2 -4) iV 2 -4 j 



(51) 



The functions X*- are integrals over the soft-gluon phase-space. In PIM kinematics one 
finds I™ = 2™' 
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(52) 



where p\ = « / 1 — 4m? /M 2 coincides with /3 in the z — > 1 limit. Furthermore x s = (1 — 
A)/(l + /3 t ), y t = -t\jm\ -1, z u = -ui/m 2 ^ - 1, and 
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lx 2 x\e 2lE 



(53) 



In 1PI kinematics one finds again that Xjf = Z^T = ^' w hile 
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In the case of 1PI kinematics, the definition of /3t should be changed to /3t = ^/ 1 — 4m? / s 
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